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We investigate the valley depolarization dynamics and valley Hall effect of exciton due to the 
electron-hole exchange interaction in mono- and bilayer M 0 S 2 by solving the kinetic spin Bloch 
equations. The effect of the exciton energy spectra by the electron-hole exchange interaction is 
explicitly considered. For the valley depolarization dynamics, in the monolayer M 0 S 2 , it is found 
that in the strong scattering regime, the conventional motional narrowing picture is no longer valid, 
and a novel valley depolarization channel is opened. For the valley Hall effect of exciton, in both the 
mono- and bilayer M 0 S 2 , with the exciton equally pumped in the K and K’ valleys, the system can 
evolve into the equilibrium state where the valley polarization is parallel to the effective magnetic 
held due to the exchange interaction. With the drift of this equilibrium state by applied uniax¬ 
ial strain, the exchange interaction can induce the momentum-dependent valley/photoluminesence 
polarization, which leads to the valley/photoluminesence Hall current. Specihcally, the disorder 
strength dependence of the valley Hall conductivity is revealed. In the strong scattering regime, the 
valley Hall conductivity decreases with the increase of the disorder strength; whereas in the weak 
scattering regime, it saturates to a constant, which can be much larger than the one in Fermi system 
due to the absence of the Pauli blocking. 

PACS numbers: 71.70.Gm, 71.35.-y, 78.67.-n, 72.25.Dc 


I. INTRODUCTION 


In recent years, as a new candidate to realize the 
valleytronics, monolayer (ML) and bilayer (BL) transi¬ 
tion metal dichalcogenides (TMDs) have attracted much 
attention^ir— To efficiently control the valley degree of 
freedom in ML and BL TMDs, both the optical^ii"— 
and electrical^i^Sr— techniques have been explored. For 
the optical method, the chiral optical valley selection 
rule allows for the optical creation of the valley po¬ 
larization, which are mainly realized by the excitonic 
excitation<^i2p— For the electrical method, due to the 
contrast Berry curvature for the electron or hole in the 
K and K’ valleys, the valley Hall effect of electron or 
hole has been predicted^ and then confirmed by the ex¬ 
periments in Mli^i and BL3^ M 0 S 2 . Furthermore, the 
method combining both the optical and electrical tech¬ 
niques to realize the valley Hall effect of trion is proposed 
theoretically.^^ This proposal is based on the fact that the 
four configurations of the trions in ML TMDs can obtain 
nonzero Berry curvature due to the electron-electron, 
hole-hole and electron-hole (e-h) exchange interactions!^ 
It can be seen that and BLi^“— TMDs pro¬ 

vide an ideal platform to study the rich valley dynamics 
based on the valley polarization or valley current. Ac¬ 
cordingly, on one hand, it is important to study the life¬ 
time of the valley polarization, i.e., the valley depolariza¬ 
tion dynamics; on the other hand, it is useful to explore 
efficient methods to create and control the valley current. 

The valley depolarization dynamics in ML and BL 
TMDs has been extensively studied, showing rich fea¬ 


tures for different membersIn ML TMDs, it has 
been experimentally found that the steady-state valley 
polarization can be large (around 50%) in M 0 S 2 , WS 2 
and WSe 2 )^'i^‘^ whereas extremely small (around 5%) 
in MoSe 2 i ^^'^^ Fast valley depolarization with the life¬ 
time about picoseconds due to the inter-valley exciton 
transition is observed It is theoretically shown that 
due to the strong Coulomb interaction, the e-h exchange 
interaction2£“— can provide an efficient valley depolar¬ 
ization channel based on the Maialle-Silva-Sham (MSS) 
mechanism j ^^ 1 ^^ Including the e-h Coulomb exchange in¬ 
teraction, the Hamiltonian of the exciton expressed by 
the center-of-mass momentum k is written as2£ 
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in which the first and second terms represent the 
kinetic energy and e-h exchange interaction, respec¬ 
tively. rriex = trie + trih is the exciton mass with me 

and rrih being the electron and hole masses; Q{k) = 
2 

2£oAt( | k | -b Ksc) 1^1^ and k± = ±iky. Here, 

Eq and K stand for the vacuum permittivity and rela¬ 
tive dielectric constant; k^c is the screening wavevector; 
(/'^°(r) = represent the exciton ground- 

state wavefunction with r being the relative coordinate 
of the electron and hole, and as denoting the exciton 
radius; oml is the material parameter.— 

This exchange interaction can cause the “precession” 
of the exciton states with a k-dependent frequency 
D(k), which causes the inhomogeneous broadening 
In analogy to the D’yakonov-Perel’ (DP) mechanismj^^ 
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in the strong scattering regime with |r2(k)|Tfe <C 1 with 
Tfe being the momentum relaxation time, the valley de¬ 
polarization time is estimated to be = (n^(k))rfc. 
Here (• • •) denotes the ensemble average. Accordingly, 
it seems that the valley depolarization should be always 
suppressed by the momentum scattering in the strong 
scattering regime. Nevertheless, Yu et al. showed that 
with the e-h exchange interaction, the energy spectra of 
the exciton is modified to be the Dirac cone<^ With this 
large modification of the exciton spectra, the momentum 
scattering should also be markedly influenced by the ex¬ 
change interaction. Although it is then pointed out that 
the Dirac cone does not exit due to the existence of the 
intra-valley e-h exchange interaction,i2S“— it is demon¬ 
strated that the exchange interaction modifies the energy 
spectra of the exciton markedlyi^^— This motivates us to 
study the valley dynamics with the exchange interaction 
explicitly modifying the energy spectra. 

Similar to the ML situation, in BL TMDs, the ex¬ 
change interaction between the four degenerate states la¬ 
beled by the valley and layer indices is also expected to 
cause the photoluminesence (PL) depolarizationi^ How¬ 
ever, it is experimentally found that different from ML 
situation, the steady-state PL polarization in the BL 
WS2 and WSe2 can be much larger than the one in ML 
under the same experimental conditions.—^— Then it is 
theoretically predicted that for the BL WS2, with the 
isotropic dielectric constant, there exits a steady state 
with the PL polarization being always half of the initial 
one due to the specific form of the exchange interaction, 
indicating that the valley depolarization time can be very 
long in BL TMDs— Specifically, the exchange interac¬ 
tion Hamiltonian between the four degenerate intra-layer 
exciton states is written as 
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noting the intra-layer relative dielectric constant; 7 = 
•\/K||/k± with K_L being the inter-layer relative dielectric 
constant. The tilde labels that the parameters in BL 
situation can be different from the ones in ML. Similar 
to the ML situation, this exchange interaction can also 
markedly modify the energy spectra of the exciton, which 
is also expected to influence the valley dynamics. More¬ 
over, the study for the PL depolarization dynamics for 
the anisotropic dielectric constant— is still lacking. 

From above analysis, it can be seen that by treating 
“valley” as “spin” in ML TMDs, the exchange interaction 
actually plays the role of the spin-orbit coupling (SOC) 
in the electronic system. In the electronic system, the 
intrinsic spin Hall effect has been well understood in the 
system with the SOC in the weak scattering limit, which 
is absent for the Rashba but can exist for other type of 
the SOC— Moreover, it is found that the intrinsic 


spin Hall conductivity is a constant in the clean sam¬ 
ple, which is nevertheless less studied in the strong scat¬ 
tering regime.— Then it is natural to expect that in 
analogy to the intrinsic spin Hall effect of electrons 
there exists the “valley” Hall effect of exciton due to the 
exchange interaction in ML TMDs. Accordingly, with 
the generation of the exciton current, which can be real¬ 
ized by applying the uniaxial strain^ the valley current 
perpendicular to the exciton current can emerge in ML 
TMDs. However, in the BL TMDs, the exchange inter¬ 
action exists between four rather than two degenerate 
exciton states, which is very different from the electronic 
system with two degenerate spin bands. It is an interest¬ 
ing problem to study whether there exists the valley Hall 
effect of exciton for the four-state system. It is empha¬ 
sized that in the previous works, the “spin” Hall effect of 
exciton has been proposed.—!^ However, it is different 
from the proposal here. In the work of Wang et al.^ the 
“spin” Hall effect arises due to the different strength of 
the SOC experienced by the electron and hole in the ex¬ 
citon; whereas in the work of Yao et al.^ it arises from 
the Berry curvature, which is in analogy to the intrinsic 
anomalous Hall effect of electron— 

In the present work, by explicitly considering exciton 
energy spectra modified by the e-h exchange interac¬ 
tion, we investigate the valley depolarization dynamics 
and valley Hall effect of exciton in ML and BL M0S2 by 
solving the kinetic spin Bloch equations (KSBEs)— For 
the valley depolarization dynamics, in the ML M0S2, it 
is found that with the exchange-interaction-modified en¬ 
ergy spectra, in the strong scattering regime, the conven¬ 
tional relation Tg oc is no longer valid. It is shown 
that a novel valley depolarization channel is opened in the 
strong scattering regime, where the valley lifetime first 
decreases and then increases with the increase of the dis¬ 
order strength, showing the Elliott- Yafel—i^ (EY) like 
behavior in the DP mechanism from the point of view of 
the spin relaxation— This channel comes from the 
inhomogeneous broadening from the module of the mo¬ 
mentum of the exciton in the exciton-disorder scattering, 
in which the same energy corresponds to different mo¬ 
mentum module with the exchange-interaction-modified 
energy spectra. This is very different from the conven¬ 
tional situation, in which the inhomogeneous broadening 
comes from the angular anisotropy of the momentum in 
the exciton-disorder scattering— * 2 ^ Moreover, due to the 
enhancement of the module-dependent inhomogeneous 
broadening by the momentum scattering, it shows EY- 
like behavior in the MSS mechanism. For the BL M0S2, 
the PL depolarization dynamics with both the isotropic 
and anisotropic dielectric constants is investigated. With 
the isotropic dielectric constant, it is shown that with 
the exchange interaction modifying the energy spectra, 
the steady state revealed in our previous worb^ still ex¬ 
ists. Whereas with the anisotropic dielectric constant, 
the steady state vanishes. However, it is found that when 
the dielectric constant is close to the isotropic situation, 
the PL polarization first decreases fast and then slowly, 
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indicating that the effective depolarization time can also 
be much longer than the ML situation. 

For the valley Hall effect of exciton, the valley Hall con¬ 
ductivity for the ML and BL M0S2 in both the weak and 
strong scattering regimes are calculated by the KSBEs. 
It is shown that with the exciton in the K and K’ valleys 
equally pumped, the system evolves into the equilibrium 
state where the valley polarization is parallel to the effec¬ 
tive magnetic field due to the exchange interaction. With 
the drift of this equilibrium state due to the applied uni¬ 
axial strain, this parallelism is broken and hence the effec¬ 
tive magnetic field can induce the momentum-dependent 
out-of-plane valley/PL polarization, which accounts for 
the valley/PL current of exciton. Furthermore, the disor¬ 
der strength dependence of the valley Hall conductivity 
is revealed. For both the ML and BL M0S2, in the strong 
scattering regime, the valley Hall conductivity decreases 
(oc t|) with the increase of the disorder strength; whereas 
in the weak scattering regime, the valley Hall conductiv¬ 
ity saturates to a constant. Specifically, it is found that 
the valley Hall conductivity in the weak scattering regime 
is proportional to the population of the exciton with zero 
momentum. By further considering that with the Bose 
distribution (therefore no the Pauli blocking), this pop¬ 
ulation can be extremely large at low temperature and 
high exciton density. Accordingly, the valley Hall con¬ 
ductivity for the exciton can be much larger than the 
one for the Fermi system. All these behaviors can be 
well understood analytically in the weak exchange inter¬ 
action limit. 

This paper is organized as follows. In Sec. m we set up 
the model and KSBEs. In Sec. imi we study the valley de¬ 
polarization dynamics and valley Hall effect of exciton in 
ML M0S2. Specifically, in Sec. lHI Al a novel valley depo¬ 
larization channel is presented; in Sec. lHIBl the disorder 
strength dependence of the valley Hall effect of exciton 
is studied first numerically and then understood analyti¬ 
cally. In Sec. m the valley depolarization dynamics and 
valley Hall effect of exciton are further discussed in BL 
M0S2. We summarize in Sec. El 


II. MODEL AND KSBES 


We start the investigation from the set up of the ki¬ 
netic equation for the exciton by considering the exciton- 
disorder scattering. We first present the effective Hamil¬ 
tonian for the exciton-disorder interaction expressed by 
the center-of-mass coordinate of the exciton, which is de¬ 
rived in Refs. [US The Hamiltonian in the disordered 
system is written as 
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where IFe(re) and Wh{rh) denote the intra-valley disor¬ 
der potential for the electron and hole, respectively, and t] 


labels the exciton state including the valley and layer in¬ 
dices. For ML M0S2, we do not consider the inter-valley 
scattering, which is suppressed for the hole with large 
splitting of the valence band, because the spin-flip scat¬ 
tering is forbidden unless the mirror reflection symmetry 
is brokeni^i^^— For BL M0S2, the inter-layer scattering 
is further neglected because only the hole with the same 
spin can hop between different layers, which is neverthe¬ 
less very weak due to the large splitting of the valen- 
cen bands, and the inter-layer hopping for the electron 
is forbidden due to the lattice symmetryJ^^— When 
the disorder is not very strong, which does not influence 
the relative motion of the exciton, the exciton-disorder 
interaction can be treated perturbatively and expressed 
by the center-of-mass coordinatei^^iffi By focusing on the 
ground state (Is-state), the center-of-mass part of the 
Hamiltonian reads 
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where R = {meVe -\- mhrh)/mex. Here, with me ~ irih in 
ML and BL M0S2, 
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describes the effective exciton-disorder interaction. From 
Eq. dS]), one notices that the charged impurity is ineffi¬ 
cient for the exciton-disorder interaction. 

Furthermore, with the e-h exchange interaction Hamil¬ 
tonian iLex(k) included in Eq. (| 3 ]), the exciton dynamics 
under the uniaxial strain can be described by the KSBEs 
including the coherent, drift and scattering terms 

^tPV. — ^tPklcoh T ^tPkldrift T '^iPklscat- (6) 

In these equations, pk represent the n x n density ma¬ 
trices of exciton with the center-of-mass momentum k at 
time t, in which the diagonal terms describe the exciton 
distribution functions and off-diagonal terms represent 
the inter-state coherence. Specifically, n = 2 and 4 for 
the ML and BL M0S2, respectively. 

In the collinear space, the coherent term is given by 

StPkIcoh =-(V^)[-f^ex(k),pk], ( 7 ) 

where [ , ] stands for the commutator. The drift term is 
denoted as 


9tPk|drift = —{F/fi) ■ VkPk, (8) 


where F represents the external force field due to the 
applied uniaxial strain. Einally, the scattering term 
i 9 tp(k, t)|scat due to the exciton-disorder scattering is 
written as 


^iPkjscat — . ^ ^ |Hk—k^l Ek,r 72 ) 


k' 7 )ir ;2 


(7k',r)i7k,r)2Pk — 7k,r;27k',r7iPk') + H-C. . (9) 
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In Eq. ([ 9 ]) 

|C/qP = J y’drdr'([C/(r)-C/o][C/(r')-C/o])e-*‘i ('-'') 

= J J drdr'Cir - = Cq, (10) 

with Uq being the average value of the disorder potential. 
C (q) is taken to be the Gaussion correlation functio n^^i^° 

C(q) = 7ry|cr|jexp(-cr|9V4), 

where Vr is the potential amplitude and an denotes the 
radius of the correlation length of the disorder. Specif¬ 
ically, when anq <C 1 , Eq. (HU actually describes the 
short-range exciton-disorder interaction. and T^.ri 

are the energy spectra of the exciton and the projec¬ 
tion matrix, whose expressions are given explicitly in Ap¬ 
pendix for both the ML and BL situations. 


III. MONOLAYER MOS2 

In this section, we investigate the valley depolarization 
dynamics and valley Hall effect for the A-exciton in the 
ML M0S2. All parameters including the band structure 
and material parameters used in our computation are 
listed in Table ID 


TABLE I: Parameters used in the computation for ML M 0 S 2 . 


me/mo 

0.35“ 

T(K) 

20 

mu/mo 

0.44“ 

Uex (cm“2) 

IQii 

K 

3.43“ 

Ue (cm~2) 

1.6 X 10 

as(nm) 

3.0 

(Tfl (run) 

6 

oml (A ) 

4.91*’ 

Vro (meV) 

4.4 X 10 

F (eV/cm) 

10 




“Reference 
Reference (2f 


In Table HI T denotes the temperature, which is 
low enough that the exciton-phonon interaction is 
neglected^^— Uex and Ue are the exciton and electron 
densities, with the former much larger than the latter. 

With these parameters, k,, = ^ ~ 

1.5 X 10 ®/m, with /ie being the chemical potential of elec¬ 
tron. 


A. Novel Valley Depolarization Dynamics 

In this subsection, we investigate the valley depolar¬ 
ization dynamics, especially focus on a novel valley de¬ 
polarization channel in ML M0S2. The valley depo¬ 
larization time is obtained by solving the KSBEs from 
the temporal evolution of the valley polarization P{t) = 


Z]k®k(0/^ex = X]kTr[pk(i)d2]/nex, with a^ being the 
^-component of the Pauli matrix. According to the chi¬ 
ral optical valley selection rulej^i^ by using the elliptically 
polarized light, the system is initialized to be 

= + (11) 

In Eq. dni), Hka = {exp[(ek - fJ.cr)/{kBT)] - 1 }“^ is 
the Bose-Einstein distribution function at temperature 
T, with Ek = ^^^^/( 2 mex) and /r-i-.j, standing for the 
chemical potentials determined by the exciton density 
?^ex=X]k 4 ’^[Pk] and the initial valley polarization P( 0 ). 
P( 0 ) = 10 % in our numerical calculation. 


1. Analytical Analysis on the Conventional Situation 


For comparison with the novel valley depolarization 
channel addressed in the next subsection fSec lIII A^ . we 
first present the analytical analysis of the conventional 
MSS mechanism^^— which has been used to under¬ 
stand the recent experimental resultsi ^^d^ As we know, 
when the splitting energy due to the exchange interac¬ 
tion, which is referred to as the “exchange energy” in 
this work, is much smaller Xhan the kinetic energy the 
exchange energy can be neglected in the energy spectra, 
i.e., Pk.r; ~ £k in Eq. Based on this approximation, 
Eq. ([ 9 ]) becomes 

27]- 

^i/^k|scat ~ E |Hk-k'P(5(ek' — ek)(pk — Pk')- (12) 

k' 


To find the valley depolarization time, we transform the 
KSBEs from the collinear space to the helix one by the 
unitary transformation pk = PkPkPk-“ In the helix rep¬ 
resentation, the KSBEs become 


dtPk + i{ef^/fi}[az,Pk] + ( 27 r/?i)a^ ^ 5 (ek' - ^k) 

k' 

X (pk — -S'kk'Pk' 5 'k'k) = 0 , ( 13 ) 

where the exchange energy = Q{k)k^ and S'kk' = 
Pj^Pk'- Here, for simplicity, we consider the situation 
with fTfl|k — k'l -C 1, hence |Pk^k'P is replaced by the 
constant = 7rHjcr|j in Eq. (fTHl) . After Fourier analysis 
with pk = one finds that the zeroth order of 

the density matrix p^ forms a closed equation, 

^ P°] + ^ - ^d^pla^ = 0 , ( 14 ) 

which is not obvious in the collinear representation!^ 
Here, 1 /r = /k? is the momentum scattering rate. 

By further noticing that S£ = Sf, one obtains from 
Eq. (HU) that 


,«(]_ 
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where Ofc = s’^lh is the precession frequency between 
different exciton states. From Eq. m, it is obtained 
that in the strong scattering limit with <C 1, 

Sk (0 ~ ^( 0 ) exp(- 4 f 2 ^rt), ( 16 ) 

and hence the valley depolarization time ~ [ 4 (n^)T]“^ 
is inversely proportional to the momentum scattering 
time, which is the motional narrowing effect in the 
random-walk theoryi^i^^— Whereas in the weak scat¬ 
tering limit with ^ 1, 

Slit) Ri P(0)e-*/(2r) cos( 2 flfct). ( 17 ) 

Hence, two factors influence the valley depolarization in 
the weak scattering limit. On one hand, the momen¬ 
tum scattering opens a valley depolarization channel due 
to the factor on the other hand, the factor 

cos(2r2fc<) can cause free induction decay due to different 
precession frequency with different momentum (inhomo¬ 
geneous broadening) 

However, with the same initial state [Eq. dill) ], above 
conventional picture obtained by the weak exchange en¬ 
ergy approximation is no longer valid when the exchange 
energy is comparable to or even larger than the kinetic 
one )^^’^^~ — which is shown in the next subsection. 

2. Momentum Scattering Dependence of the Novel Valley 
Depolarization 

In ML M0S2, due to the strong Coulomb interaction 
and the large exciton mass, the exchange energy is com¬ 
parable to the kinetic onei ^^i^^~ — This is true even when 
the screening effect due to the residue electron in the 
sample is considered (Table H]). Therefore, the exchange 
energy should enter the energy spectra in the scatter¬ 
ing term when calculating the valley depolarization time. 
Consequently, a novel valley depolarization channel in 
the strong scattering regime, in which the valley depo¬ 
larization is enhanced rather than suppressed by the mo¬ 
mentum scattering, showing the EY-like behavior^^^i^ is 
switched on. 

To understand the new valley depolarization channel, 
we focus on a simplified model, where only the diagonal 
elements in the projection matrix [Eq. dA 2 l) ] are retained. 
The corresponding scattering term reads 

^tPklscat ~ ^ ^ / d£/c'd(ii/k', 7 )i ^'k,r/2)(pk Pk'\ 

( 18 ) 

where P% = ^ ^ f^^k'Pk'- Furthermore, by means of 
the Fourier analysis, the density matrix is expanded as 
Pk = Pfc + 'l 2 i^o where the zeroth and non-zeroth 

components are written separately. It can be seen that 
for the zeroth order of pk, the scattering term 

l 9 tPk|scat ~ ^ ' j d£k'^{Eu',ui ~ Eu,r]2)iPk ~ Pk') 

(19) 


is nonzero for the exciton-disorder scattering between dif¬ 
ferent energy branches; whereas it is forbidden in the 
conventional situation from Eq. uni). 

Therefore, Eq. m opens an additional valley depo¬ 
larization channel by causing the module-dependent in¬ 
homogeneous broadening. This is very different from 
the conventional situation, in which the inhomogeneous 
broadening arises from the angular anisotropy of the mo¬ 
mentum in the exciton-disorder scattering^i*^ With this 
enhancement of the inhomogeneous broadening, the val¬ 
ley depolarization tends to be enhanced. However, this 
additional channel also enhances the momentum scatter¬ 
ing, which tends to suppress the valley depolarization. 
Therefore, there exists the competition between the ef¬ 
fective inhomogeneous broadening and momentum scat¬ 
tering in this new channel. It is demonstrated that the 
EY-like behavior exactly comes from this enhancement 
of the inhomogeneous broadening. It is shown in Eig. [ 1 ] 
that compared to the full calcualtion (the red solid curve 
with circles), when the additional valley depolarization 
channel is removed, the EY-like behavior vanishes in the 
pink chain curve. 



FIG. 1: (Color online) Disorder strength dependence of the 
valley depolarization time in ML M 0 S 2 . The boundary be¬ 
tween the weak and strong scattering regimes [2(flfc)r « 1] 
is denoted by the vertical cyan dashed line. The valley de¬ 
polarization times with both the short-range scattering (the 
red solid curve with circles) and Gaussion correlation func¬ 
tion with aR = 6 nm (the blue dashed curve with squares) 
are shown. With short-range scattering, the red solid curve 
with circles and gray dashed curve correspond to the valley 
depolarization times fitted from 5Dk and |Sk|, respec¬ 
tively. For comparison, the conventional situation is plotted 
by the blue (JDk ®k) and orange (JDj^ |Sk|) dashed curve with 
crosses. Finally, the green chain curve corresponds to the sim¬ 
plified model [Eq. (1181) 1. and the pink chain curve represents 
the situation with Eq. (O removed from Eq. ®- 

Figure [T] shows the disorder strength dependence of the 
valley depolarization time computed, based on the ma- 
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terial parameters shown in Table HI In Fig. [ 1 ] the bound¬ 
ary between the weak and strong scattering regimes 
[iyuiymY ~ 0 . 06 ] is shown as the vertical cyan dashed 
line. Accordingly, it is shown that with both the short- 
range (the red solid curve with circles) and Gaussian cor¬ 
relation funciton (the blue dashed curve with squares) in 
the scattering term, in the weak scattering regime, the 
valley depolarization time increases with the increase of 
the disorder strength; whereas in the strong scattering 
regime, the valley depolarization is first enhanced and 
then suppressed by the momentum scattering, which is 
referred to as the EY-like and normal strong scattering 
regimes. Furthermore, these two curves for the short- 
range and Gaussion correlation function coincide with 
each other, showing that the short-range scattering is 
a good approximation for the exiton-disorder scattering 
here. This is because for the Bose-Einstein distribution 
at low temperature here, qan <C 1 is satisfied, and hence 
the exponential function in Eq. dm) can be neglected. 

For comparison, the conventional situation is also com¬ 
puted, shown by the blue dashed curve with crosses in 
Fig. [TJ It can be seen that with the increase of the dis¬ 
order strength, the valley depolarization time increases 
monotonically in both the weak and strong scattering 
regimes. In the strong scattering regime, this confirms 
Eq. (fTBl) where the valley depolarization is suppressed by 
the momentum scattering. Whereas in the weak scat¬ 
tering regime, from Eq. (HU, it seems that the valley 
depolarization should be enhanced by the momentum 
scattering. However, because the energy dispersion in 
the Bose-Einstein distribution, the free induction decay 
dominates the valley depolarization, which can be sup¬ 
pressed by the momentum scattering. To see this point, 
we plot the valley depolarization time r* fitted from in¬ 
coherently summed spin polarization | Sk | where 
the free induction decay is destroyed. It is shown by the 
orange dashed curve with crosses that with the increase 
of the disorder strength, r* decreases in the weak scatter¬ 
ing regime because the momentum scattering can directly 
open a channel for the valley depolarization in 

Eq. (flTl) ]. but increases in the strong scattering regime. 

It is interesting to see that the behavior of the valley 
depolarization with the exchange-interaction-modified 
energy spectra is similar to the conventional situation 
in the weak scattering regime, but very different in the 
strong scattering regime. In the weak scattering regime, 
the new valley polarization channel is not important, 
which can be seen from Fig.[I]that the valley depolariza¬ 
tion with (green chain curve) and without (pink chain 
curve) this channel almost coincides with each other. 
Therefore, same as the conventional situation, due to the 
suppression of the free induction decay by the scattering, 
the valley depolarization time increases with the increase 
of the disorder strength. By further fitting r* from inco¬ 
herently summed spin polarization |Sk| (gray dashed 
curve), one observes r* decreases with the increase of 
the disorder strength in the weak scattering regime due 
to the destroy of the free induction decay. In the EY-like 


regime, the enhancement of the valley depolarization by 
the scattering comes from the enhancement of inhomoge¬ 
neous broadening due to the novel valley depolarization 
channel [Eq. (fTOl) ]. Finally, in the normal strong scatter¬ 
ing regime, the momentum scattering is very strong. The 
enhancement of the momentum scattering in Eq. m 
becomes more important than the enhancement of the 
inhomogeneous broadening, and hence the valley depo¬ 
larization is suppressed by the momentum scattering. 

Finally, it can be seen from Fig.[I]that the green chain 
curve calculated from the simplified model [Eq. (fT51) ] al¬ 
most coincides with the one by full calculation (the red 
solid curve with circles). Therefore, it seems that the 
off-diagonal elements in the projection matrix play less 
important role in the valley depolarization in the regimes 
we study. However, it influences the behavior of the tem¬ 
poral evolution of the valley polarization in the normal 
strong scattering regime, leading to the oscillations of the 
valley polarization (refer to Appendix |B]) . 

B. Valley Hall Effect of Exciton 

In this part, we study the valley Hall effect of exciton 
in ML M0S2 both numerically and analytically. In the 
calculation, the initial state is set to be the equally popu¬ 
lated Bose-Einstein distribution in the K and K’ valleys, 
i.e., 

Pk(0) = B°i, (20) 

where is the Bose-Einstein distribution function. The 
computation parameters are list in Table [H With this 
experimentally-realized initial state^i^ we first show that 
the system can evolve to the equilibrium state in which 
the “spin” vectors [S^ = Tj^ipk^x), = Tr(pkdy) and 
S^] are parallel to the k-dependent magnetic field due 
to the exchange interaction [Sec. IHI B Ij . Then we show 
that with this equilibrium state, after applying the exter¬ 
nal force field due to the uniaxial strain^i^ the drift of 
this equilibrium state can induce the valley Hall current 
of exciton ISec. IHI 5 ^ . 

1. Equilibrium State without External Eorce Eield 

Before the concrete study of the valley Hall effect of 
exciton in ML M0S2, it is important to know the prop¬ 
erty of the equilibrium state with the exciton initially 
equally-populated in the K and K’ valleys. In Figs, ^a.) 
and (b), it is shown that after a long time, the sys¬ 
tem evolves into the equilibrium state which corresponds 
to the “spin” separation for and in the momen¬ 
tum space, respectively. Specifically, in this equilibrium 
state, the spin vectors and are parallel to the k- 
dependent magnetic field due to the exchange interaction 
= — 2 e“cos( 20 k) and = — 2 e|’^sin( 20 k)- Below 
we analytically demonstrate this property in the weak 
exchange interaction limit. 
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FIG. 2: (Color online) Momentum distribution of the “spin” 
vectors (a) and (b) in the equilibrium state, fco « 
7.9 X 10^/m is the Fermi wave-vector of the system. The cal¬ 
culation shows that in the equilibrium state, oc cos(20k) 
and oc sin(26'k), which are parallel to the in-plane ef¬ 
fective magnetic field due to the exchange interaction along 
i-direction [12^ = —2e^* cos(26k)] and y-direction [12^ = 
—2£™ sin(26'k)], respectively. 


In the weak exchange interaction limit, we expand the 
energy spectra in the linear order of the exchange en¬ 
ergy in the scattering term [Eq. (jH])], and then derive 
the equilibrium state to be (the derivation is referred to 
Appendix ICl) 

pi^Bli + Jf^^dBt/dsk. ( 21 ) 

Here, with the diagonal elements contributing to the en¬ 
ergy spectra, the exchange interaction Hamiltonian 

f 0 —\ 

\ j ( 22 ) 

only contains the off-diagonal elements. This equilibrium 
state corresponds to the spin vectors 


Obviously, the spin vectors and 5"^ are parallel to the 
k-dependent magnetic field 12^ and j2^, respectively. 


2. Valley Hall Effect of Exciton 


In this subsection, we study the valley Hall effect of 
exciton. Specifically, we explicitly show the disorder 
strength (momentum scattering) dependence of the val¬ 
ley Hall conductivity From the initial state Eq. (EOI), 
we numerically calculate the steady state density matrix 
with the applied field after long temporal evolution by 
the KSBEs. Then the valley Hall conductivity is calcu¬ 
lated. With the valley Hall current defined as 


fx = + VyU^)\^ = cr^E/|e|, (25) 


the valley Hall conductivity is expressed as 


= 


2F 


Tr p\^{a^Vy + Vy&z) 


(26) 


Here, |e| is the electron charge and Vy = hky/niex + 
dH^^{'k)/dky is the velocity operator. 

The results are shown in Fig. [31 With the weak ex¬ 
ternal force field, the system is in the linear regime, as 
shown in Fig. [3] that the momentum scattering time re¬ 
veals linear dependence on the disorder strength. In this 
regime, it can be seen from Fig.|3]that in the strong scat¬ 
tering regime, the valley Hall conductivity decreases with 
the increase of the disorder strength, showing the depen¬ 
dence oc T^; whereas in the weak scattering regime, 
the valley Hall conductivity saturates to a constant. Be¬ 
low we show analytically that these features in the mo¬ 
mentum scattering dependence can be well understood 
in the weak exchange interaction approximation. 

Here, we outline the main results to obtain the physi¬ 
cal picture of the valley Hall effect of exciton. It has been 
shown that without the external force field, the density 
matrix in the equilibrium state is written as Eq. m, 
which commutes with the exchange interaction Hamilto¬ 
nian. After applying the external force field, the density 
matrix in the steady state is derived in Appendix |D1 as 
shown by Eq. (ID5I) . It is shown that based on the equi¬ 
librium state [Eq. (|2T]) ]. after applying the external field, 
the drift part of the density matrix is [Eq. (ID2I) ] 



Ft d 
h dkx 




dsk 


{27) 


Obviously, this drift density matrix no longer commutes 
with the exchange interaction Hamiltonian, i.e., it can in¬ 
duce the momentum-dependent “spin” polarization along 
the z-direction (valley polarization). This can be seen as 
follows. In Eq. (iDSl) . the induced density matrix respon¬ 
sible for the valley polarization is written as 


(23) 

(24) 


Pk 


n? dsk 1 + 402r2 ^ ^ ^ dkx 


(28) 


= -2ercos(2dk)aH°/9efe, 
SI = -2eTsin{2e^)dBlldek. 
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FIG. 3: (Color online) Disorder strength dependence of the 
valley Hall conductivity and the momentum scattering time 
(note the scale is on the right hand side of the frame). The 
cyan dashed line labels the boundary between the weak and 
strong scattering regimes. The green chain (short-range) 
and orange dashed (Gaussian correlation) curves with squares 
show that the momentum scattering time is linearly depen¬ 
dent on the disorder strength. For the valley Hall conduc¬ 
tivity, it can be seen from the red solid curve with squares 
(short-range) and blue dashed curve with circles (Gaussian 
correlation) that in the strong scattering limit, the valley 
Hall conductivity decreases with the increase of the disorder 
strength with the dependence oc whereas in the weak 
scattering regime, the valley Hall conductivity saturates to a 
constant. 


Then with the exchange interaction Hamiltonian 
[Eq. (|2^ ]. it can be obtained that 


~ + 2Qkyay-, (29) 

K. MQ^k^kyU ^. (30) 


In the derivation, we have used the fact that when the 
screening effect due to the residue electron is considered, 
Q(/c) is approximately a constant with our computation 
parameters (Table |T|. From Eq. (1501) . one observes that 
the induced density matrix is proportional to and de¬ 
pendent on the momentum ky. 

Then the valley Hall conductivity can be calculated. 
With the definition of the valley Hall conductivity 
[Eq. (HH)], only the third term on the right-hand side 
of Eq. (ID5I) . i.e., Eq. (l28l) . contributes to the valley Hall 
conductivity. From Eq. dUl), one has 


jej p ^ dBl 

~ hi + 


(31) 


From above equation, one finds that in the strong scat- 
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z 

X 


h 



dB'l 2 2 
dsk-^flkT^, 
d£k 


(32) 


which is proportional to r^. Whereas in the weak scat¬ 
tering regime with HfcT > 1, 


^ MhLo, (33) 

which is independent on the exchange interaction 
strength and momentum scattering. Specifically, one 
observes that is proportional to which can 

be extremely large when the system is close to the 
Bose-Einstein condensation. Here, with the computa¬ 
tion parameters (Table |T|, ~ hence ~ 

0 .6|e|//i, which gives a good estimate to the calculated 
one 0A\e\/h with strong exchange interaction. Moreover, 
from Eq. ((551) . one observes that with the higher exci- 
ton density and/or lower temperature, is large and 

hence the valley Hall conductivity. According to our cal¬ 
culation, with the exciton density riex = 5 x 10^^ cm“^ 
at 10 K, ~ 5.5|e|//i in the weak scattering limit. This 
is much larger than the one in Fermi system, with the 
latter being limited by the Pauli blockingj ^^d^ 

Finally, we summarize the physical picture of the val¬ 
ley Hall effect of exciton as follows. First of all, it is 
understood that in the equilibrium state, the “spin” vec¬ 
tor of any momentum k is parallel to the k-dependent 
effective magnetic field due to the e-h exchange interac¬ 
tion (Sec. IHIB ip . Then by applying the force field due 
to the uniaxial strain)^ the “spin” vector is no longer 
parallel to the effective magnetic field. Accordingly, the 
“spin” vector can rotate around the effective magnetic 
field, and the momentum-dependent valley polarization 
is induced. Specifically, for the exciton with opposite ky, 
the ai-component of the effective magnetic field is along 
the opposite direction, and hence the induced valley po¬ 
larization is also opposite. Consequently, the valley cur¬ 
rent perpendicular to the driven exciton current is estab¬ 
lished. 

It is emphasized that the physical picture addressed 
above is in analogy to the intrinsic spin Hall effect of the 
electron i^®d^~ — Nevertheless, two new features in this 
mechanism are further revealed here. On one hand, it 
is revealed that in the dirty sample corresponding to the 
strong scattering regime, the intrinsic “spin” Hall effect is 
markedly suppressed by the momentum scattering, with 
its conductivity proportional to r^. On the other hand, 
in the weak scattering regime, the Bose system with no 
Pauli blocking provides an ideal platform to realize large 
“spin” Hall conductivity, which can be much larger than 
the one in the Fermi system, especially when the system 
is close to the Bose-Einstein condensation. 









9 


IV. BILAYER MOS 2 

In this section, we investigate the valley depolariza¬ 
tion dynamics and valley Hall effect for the A-exciton, 
which is four-fold-degenerate, in BL M 0 S 2 . In our pre¬ 
vious work, a steady state in the PL depolarization dy¬ 
namics in BL WS 2 with the isotropic dielectric constant 
was revealed in the situation without the energy spectra 
modified by the exchange interaction.^^ However, with 
the exchange-interaction-modified energy spectra, as re¬ 
vealed in ML M 0 S 2 (Sec. IBII) . the valley dynamics be¬ 
comes very different from the conventional situation and 
the valley Hall effect of exciton can arise. So far, the PL 
depolarization dynamics for BL TMDs with anisotropic 
dielectric constant is still lacking. These motivate us 
to calculate the related PL dynamics for the four-fold- 
degenerate states in BL system with the energy spectra 
modified by the exchange interaction. All parameters in¬ 
cluding the band structure and material parameters used 
in our computation are listed in Table m 


TABLE II: Parameters used in the computation for BL M 0 S 2 . 


me/mo 

0.35“ 

r(K) 

20 

mh/mo 

0.44“ 

Sex (cm”'^) 

10“ 

K\\ 

4.8'’ 

Ue (cm”'^) 

3.3 X 10 “ 

as (nm) 

3.0 

an (nm) 

6 

obl (A ) 

4.51“ 

Vro (meV) 

5.4 X 10"^ 

F (eV/cm) 

10 




a 

b 

c 


Reference 

Reference 

Reference 


6E| 

^1 

34| 


It is emphasized that in BL TMDs, the centre-inversion 
symmetry exists Hence, no valley polarization 
but the PL polarization can be created by the chiral opti¬ 
cal valley selection ruleJ^^— Similarly, no valley current 
but the PL current can be created by the valley Hall ef¬ 
fect due to the symmetry, which can be measured at the 
edges of the device channel in the experiment 


A. PL Depolarization Dynamics 

In this part, we focus on the new feature of the PL 
depolarization dynamics arising in BL system compared 
to the ML situation. The PL depolarization dynamics is 
obtained by solving the KSBEs from the temporal evo¬ 
lution of the PL polarization P{t) = 
with 


(1 

0 

0 

o \ 

0 

-1 

0 

0 

0 

0 

-1 

0 


0 

0 

1 / 


(34) 


In the calculation, the initial condition is set to be 

-Bkt + Bk4. Bkt - Bkj. - , . 

/3k(0) =-^ -b-- Iz. (35) 

P{0) is set to be 10% in our numerical calculation. Be¬ 
low, the PL depolarization dynamics with isotropic and 
anisotropic dielectric constants are investigated, respec¬ 
tively. 

For the isotropic dielectric constant (7 = 1), it is 
shown in Fig. |4] that no matter the system lies in the 
weak or strong scattering regime, there is always a steady 
state with the PL polarization being half of the ini¬ 
tial one. This is the same as our previous prediction 
without the exchange-interaction-modified energy spec¬ 
tra (the red solid curve in Fig. |3])i^ However, when the 
exchange interaction markedly modifies the energy spec¬ 
tra, the density matrix in the steady state is found to be 
different from the previous onei^ 



x(ps) 

FIG. 4: (Color online) Temporal evolution of the PL polariza¬ 
tion in the BL M 0 S 2 with isotropic dielectric constant (7 = 1 ) 
and different disorder strength. It is shown that no mat¬ 
ter the disorder strength is weak or strong, there is always a 
steady state with the PL polarization being half of the initial 
one. The red solid curve is calculated without the exchange- 
interaction-modified energy spectra.— 


In our previous work, with the initial condition 
Eq. (l35|) . when the system lies in the steady state, the 
density matrix is found to be^ 


B, 


Pk = 


kt ■ 


B, 


k^ 


Bkt — B 


kj. 


/ 1 0 0 1 \ 
0 - 1-10 
0 - 1-10 
V1 0 01/ 


■ (36) 


Here, with the exchange interaction markedly modifies 
the energy spectra, the form of the density matrix in 
the steady state is different, in which no zero elements 
arise and hence all states are correlated to each other. 
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This can be understood in the weak exchange interaction 
approximation. As a simplified model, with the diagonal 
and off-diagonal elements entering the energy spectra of 
exciton, the effective exchange interaction in BL M 0 S 2 is 
written as 


_g2*ek Q \ 

0 0 

0 0 ’ 

Q j 

(37) 

where = Q(k)k^. In the weak exchange interaction 
approximation, the KSBEs for the BL M 0 S 2 are similar 
to the ML situation (Appendices O and [D]) , 


JUy. 


^ex 


-220k 


7e 
—e 


V 0 


0 

-210^ 


^jOk 

dt 


F dp\^ 
h dkx 


+ ^[■^k^Pk] 


Pk-P° _ ^ 
r T 


X dS{ek' - Pk - Pk'} = 0, 


d0k^ 

( 2^)2 

(38) 


where {,} denotes the anti-commutator. Without the 
applied field, in the steady state, the first three terms 
in the left-hand side of Eq. (1551) are zero. Based on the 
conventional density matrix in the steady state [Eq. (1551) ] 
and following the iteration technique introduced in Ap¬ 
pendix [CJ one obtains the steady state here (7 = 1), 


Pk ~ Pk + 

= Pi + (l/2)j4«"(aBkt/aek + aBk^/ask). (39) 

Obviously, Eq. (1551) commutes with the exchange inter¬ 
action Hamiltonian Eq. ((37|) . Specifically, due to in 
Eq. (|39)) . all states become correlated to each other. 

For the anisotropic dielectric constant, it is seen from 
Fig. [5] that compared to the isotropic case 7 = 1, when 
the dielectric constant is tuned to be anisotropic with 
7 = 1.1 (the blue chain curve) and 1.2 (the green dashed 
curve), the steady state vanishes. However, when 7 is 
close to the isotropic situation, the PL polarization first 
decreases fast and then slowly. Accordingly, the effective 
depolarization time can also be much longer than the ML 
situation, shown as 7 = 0 by the orange dashed curve. 


B. Valley Hall Effect of Exciton 

In this subsection, we investigate the valley Hall effect 
of exciton in BL M 0 S 2 . In the calculation, the initial 
state is set to be the equally populated Bose-Einstein 
distribution in the K and K’ valleys in both the upper 
and lower layers, i.e.. 


Pk(0) = 5^4X4. (40) 

From the KSBEs, with the applied force field, the steady- 
state density matrix is calculated and then used to cal¬ 
culate the valley Hall conductivity. In analogy to the 
ML situation, with the PL current which carries the PL 



FIG. 5 : (Color online) Temporal evolution of the PL polar¬ 
ization in the BL M0S2 with anisotropic dielectric constant. 
When the dielectric constant is tuned to be anisotropic with 
7 = 1.1 (the blue chain curve) and 1.2 (the green dashed 
curve), the steady state vanishes. For comparison, the time 
evolution of the valley polarization in ML M0S2 is also plot¬ 
ted, which is shown as 7 = 0 by the orange dashed curve. 


polarization defined as 


jx Pk^( 4 tly + Vyl^) = CT^E/jej, ( 41 ) 


the valley Hall conductivity in BL M 0 S 2 is expressed 
as 


( 42 ) 


Here, Vy = hky/fhe-x + dF[f^{k) / dky. The calculated 
results with both the isotropic and anisotropic dielectric 
constants are summarized in Fig. |51 

In Fig. ini the disorder strength dependence of the val¬ 
ley Hall conductivity and momentum scattering time is 
plotted. The curves for the momentum scattering time 
show that the system lies in the linear regime. It can be 
seen from Fig. [5] that in this regime, no matter the di¬ 
electric constant is isotropic with 7 = 1 (the red solid 
curve with squares) or anisotropic with 7 = 1.2 (the 
blue dashed curve with squares), in the strong scattering 
regime, the valley Hall conductivity decreases with the in¬ 
crease of the disorder strength, showing the dependence 
oc T^; whereas in the weak scattering regime, the val¬ 
ley Hall conductivity saturates to a constant. This shows 
that although the PL depolarization dynamics in the BL 
M 0 S 2 is different from the ML situation, the properties 
of the valley Hall conductivity between them are simi¬ 
lar. Below, the valley Hall conductivity is also derived in 
the weak exchange interaction approximation, with the 
steady-state density matrix with the applied field derived 
in Appendix [P] 
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V. SUMMARY 
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FIG. 6 : (Color online) Disorder strength dependence of the 
valley Hall conductivity in BL M 0 S 2 . The blue chain (7 = 1) 
and orange dashed (7 = 1 . 2 ) curves with squares represent 
the momentum scattering time, showing the system lies in 
the linear regime with the applied field (note the scale is on 
the right hand side of the frame). It is shown that no matter 
the dielectric constant is isotropic with 7 = 1 (the red solid 
curve with squares) or anisotropic with 7 = 1.2 (the blue 
dashed curve with squares), in the strong scattering regime, 
the valley Hall conductivity decreases with the increase of the 
disorder strength, showing the dependence oc whereas 
in the weak scattering regime, the valley Hall conductivity 
saturates to a constant. 


As a simplified model, the exchange interaction Hamil¬ 
tonian Eq. (1571) is used. From Eq. (1551) . it can be seen 
that only the third term on the left-hand side of Eq. (ID5I) 
contributes to the valley Hall conductivity. In Eq. (IDSI) . 
with the exchange interaction Eq. (EZl), one finds 




AiQ'^k^ky 


I ( 7 ^ + 1 ) 0 

0 -(7^ + 1) 

0 27 

V -27 0 


0 

27 

-(7^ + 1) 
0 


-27 \ 

0 

0 

(7^ + 1)/ 


It is noted that here when the screening effect is consid¬ 
ered, Q{k) has been treated as a constant. Accordingly, 
the valley Hall conductivity in BL M 0 S 2 is written as 


cr 


Z 

X 



dBl 

defc 1 + ’ 


(44) 


where Vlk = Obviously, for the valley Hall conduc¬ 

tivity in BL M 0 S 2 , it is interesting to see that Eq. (H51) 
is similar to Eq. dSID in ML M 0 S 2 . Therefore, in the 
weak and strong scattering regimes, similar features for 
the valley Hall conductivity to the one in ML situation 
can be obtained, as addressed in Sec. IHIB2I 


In summary, we have investigated the valley depo¬ 
larization dynamics and valley Hall effect of exciton in 
ML and BL M 0 S 2 by solving the KSBEsi^ The effect 
of the exchange-interation-modified energy spectra is ex¬ 
plicitly considered. For the valley depolarization dynam¬ 
ics, in ML M 0 S 2 , it is interesting to find that the con¬ 
ventional motional narrowing relation Tg oc in the 
strong scattering regime is no longer valid. It is revealed 
that in this regime, a novel valley depolarization channel 
is opened, where the valley lifetime first decreases and 
then increases with the increase of the disorder strength, 
showing the EY-lik o®^’®^ behavior from the point view 
of the spin relaxation.— This channel comes from the 
newly module-dependent inhomogeneous broadening in 
the exciton-disorder scattering, in which the same en¬ 
ergy corresponds to different momentum modules due to 
the exchange-interaction-modified energy spectra. This 
is very different from the conventional situation, in which 
the inhomogeneous broadening comes from the angu¬ 
lar anisotropy of the momentum in the exciton-disorder 
scattering—!^ Moreover, due to the enhancement of the 
inhomogeneous broadening by this channel, EY-like be¬ 
havior arises in the MSS mechanism. 

For BL M 0 S 2 , the PL depolarization dynamics with 
both the isotropic and anisotropic dielectric constants is 
investigated, which are found very different from the ML 
situation. With the isotropic dielectric constant, it is 
shown that with the exchange-interaction-modified en¬ 
ergy spectra, the steady state revealed in our previous 
worb^i still exists. Whereas with the anisotropic dielec¬ 
tric constant, the steady state vanishes. However, it is 
found that when the dielectric constant is close to the 
isotropic situation, the PL polarization first decreases 
fast and then slowly, indicating that the effective depolar¬ 
ization time can be much longer than the ML situation. 

For the valley Hall effect of exciton, the valley Hall 
conductivity for ML and BL M 0 S 2 in both the weak 
and strong scattering regimes are studied numerically 
and analytically. We show that with the exciton equally 
pumped in the K and K’ valleys, the exciton states evolve 
into the equilibrium state with the valley polarization 
parallel to the momentum-dependent effective magnetic 
field due to the exchange interaction. Then with the 
drift of this equilibrium state due to the applied uniaxial 
strain, this parallelism is broken and hence the effective 
magnetic field can induce the momentum-dependent val¬ 
ley/PL polarization, which accounts for the valley/PL 
current. This mechanism is in analogy to the intrinsic 
spin Hall effect of the electron— 

Furthermore, it is found that althougth the valley/PL 
depolarization dynamics is very different between the 
ML and BL situations, the valley Hall effect shows sim¬ 
ilar features in the momentum scattering dependence. 
Specihcally, in the strong scattering regime, the valley 
Hall conductivity decreases with the increase of the dis¬ 
order strength (oc r^); whereas in the weak scattering 
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regime, the valley Hall conductivity saturates to a con¬ 
stant, which is proportional to the population of exciton 
with k = 0. Therefore, on one hand, in the dirty sam¬ 
ple corresponding to the strong scattering regime, the 
valley Hall effect is markedly suppressed by the momen¬ 
tum scattering; on the other hand, in the weak scattering 
regime, the Bose system with no Pauli blocking provides 
an ideal platform to realize large “spin” Hall conductiv¬ 
ity, which can be much larger than the one in the Fermi 
system, especially when the system is close to the Bose- 
Einstein condensation. 
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Appendix A: Energy spectra and projection matrix 


In this appendix, we present the energy spectra and 
projection matrix for ML and BL M 0 S 2 . For the ML 
situation, the energy spectra read 

^ ;i2^2/(2mex) + Q{k)k^ ± Q{k)k'^. (Al) 

The projection matrices are 


Tkl] 

2fc2 ^ ^ki e j ■ 


(A2) 


For the BL situation, the energy spectra are 

Ell = Ell = ft"fcV(2mex), 

Ell = + 2(1 - i)Q{k)k\ 

Ell = ?i"fcV(2mex) + 2(1 + j)Q{k)k‘^. (A3) 


The corresponding projection matrices are given by 


-^k,l 


rr^BL 


rpBL _ ^ 

k.3 4^2 



/ 0 0 0 0 \ 

1 

0 0 


0 0 0 0 


\0 kl 0 y 


/ /c 2 0 o\ 

1 

0 0 0 0 


kl 0 k^ 0 


0 

0 

0 

0 

/ e 

r .2 1,2 r 2 

r\j 1 ^-j- ^ 

-kl 

k^ k^ —kl 

-kl 

—kl 


)l 2 r 2 r 2 

l\i 1 ^ 


(A4) 


(A5) 


(A6) 


and 


rr-iBL _ 

-^k^ — 



to 

+10 

1 

1 

to 

1 

kl e -k^ -kl 

4F 

—kl —kl 


\ -e -kl kl e / 


(A7) 


Appendix B: Role of off-diagonal elements of 
Eq. (IA2I) on valley depolarization 

Here, we address the role of the off-diagonal elements 
in the projection matrix [Eq. (IA2I) ] on the temporal evo¬ 
lution of valley polarization in ML M 0 S 2 . It is shown 
in Fig. 13 by the dashed curves that without the off- 
diagonal elements in the projection matrix, the valley 
depolarization is a little enhanced compared to the full 
calculation by the solid curves in the weak scattering 
[{V/VroY = 0.01], EY-like [(F/Vro)^ = 0.5] and normal 
strong scattering [(F/Vro)^ = 10] regimes. Moreover, in 
the normal strong scattering regime, the off-diagonal el¬ 
ements of Eq. (IA2p cause the oscillations in the temporal 
evolution of the valley polarization, shown as the black 
solid curve. By removing the off-diagonal elements in the 
projection matrix, the oscillation vanishes and the valley 
polarization becomes just the oscillation exponential de¬ 
cay (the blue dashed curve). 



1 2 3 4 5 6 7 8 

T(ps) 

FIG. 7: (Color online) Temporal evolution of the valley polar¬ 
ization in ML M 0 S 2 with (solid curves) and without (dashed 
curves) the off-diagonal elements in Eq. (IA2I) . The role of 
the off-diagonal elements in the weak scattering [(F/Vro)^ = 
0.01], EY-like [(F/Vro)^ = 0.5] and normal strong scattering 
[{V/VroY ~ 10] regimes are shown. 


Appendix C: Analysis on the equilibrium state 

We focus on the situation where the exchange energy 
is much smaller than the kinetic one. With the weak 
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exchange interaction the delta function in Eq. ([9]) 

is expanded, e.g., 

S{E^,^+ - Ek.+ ) = (5(£k' - fk + 4' - eD 
^ ^(ek' - Sk) + - en (Cl) 

With the linear order of the exchange energy retained in 
the scattering term [Eq. ([5])], the KSBEs are written as 

^ -^k)iPk' -Pk) 

k' 

k' 

In the equilibrium state, dtp^ = 0 and [j^^, Pk] = 0. 
Hence, one obtains 

Pk = P° + / ^dSie,, - - Jrk'=,^ pl - p^,}. 

(C3) 

This integral equation can be approximately solved by 
using the iteration technique. It is assumed that Pk = 
+ J2^=i Pk"^- substituting S^d + Pk ^ on the left- 
hand side and B'^I on the right-hand side of Eq. (IC3p . one 
obtains = M’^dB'^jdEk- By repeating this process, 

one finds p^ is proportional to (e^)^. Here, we only 
keep the linear order in the exchange energy, i.e., 

pl^Bli + M^-dBlldsk. (C4) 

Obviously, p^ commutes with . 


Eq. (jPlI) is an integral-differential equation, which can 
be solved by the iteration technique approximately. The 
density matrix is assumed to be Pk = Pk + Pk"^ 

with pI^"^ oc 

The zeroth order of Eq. dnu is exactly Eq. (IC3I) . whose 
solution has been expressed by Pk [Eq. (|2T]) ]. The first 
order of Eq. m reads 




h dk^ 


^ (R0f< 


(D2) 


which is just the drift form of the equilibrium state. 


One notes that the drift density matrix p^^ no longer 
commutes with which causes the precession of the 

“spin” vectors around the k-dependent effective magnetic 
field. The n-th order in > 2) density matrix satisfies 


Fdp\ 


n-l) 


h dkx 






2^ex 
jTy. , 


(n-l) 

k 

dck 


dP] 




x{^k'^pL”■'^} = 0, 


1 r (n) -(n—1)] 

-K -Pic J 
dPk' 


27r 


■d6{ek' - Ek) 


(D3) 


where = l/(27r) f dPkPk"^- is complex, but 

fortunately it can be much simplified if only the density 
matrix in the linear order of F is retained (linear regime). 
Furthermore, p^^^ (n > 1) and the last term on the left- 
hand side of Eq. (ID3p are exactly zero due to the angle 
integration. Finally, one obtains (n > 2), 


(z/n) [^k"^ + pEEt = 0. (D4) 


Appendix D: Solution of KSBEs with an applied 
field 


When the exchange interaction is weak, in the steady 
state (i9tpk = 0), the KSBEs with the external force field 
can be simplified to be 


F dpi, 
h dkx 


+ ^ Pk] + - (pk — Pk) 


1 



} 


TT f dEk'dPi,, d5{Ek' - Ek) r ,^ex T r, 

-7/-(SF— 


(Dl) 


With p^E^ [Eq. (ID2p ] known, p^"^ (n > 2) can be ob¬ 
tained. By summing Pk”\ one comes to a dosed form 
of the density matrix for ML M 0 S 2 (for BL M 0 S 2 , one 
replaces by and ^k by Hfc), 


Pk 




9B°k\ 

15f> 


F i dBl 

h dkx k? dEk 1 -I- 4H^t^ 


X 


[■^k 


d^ 

dkx 


F djr^^ dBl T 
h dkx dEk 1 + 


(D5) 
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